BCS-BEC crossover at finite temperature for superfiuid trapped Fermi atoms 
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We consider the BCS-BEC crossover for a system of trapped Fermi atoms at finite temperature, both 
below and above the superfiuid critical temperature, by including fiuctuations beyond mean field. 
We determine the superfiuid critical temperature and the pair-breaking temperature as functions of 
the attractive interaction between Fermi atoms, from the weak- to the strong-coupling limit (where 
bosonic molecules form as bound- fermion pairs). Density profiles in the trap are also obtained for 
all temperatures and couplings. 



PACS number(s): 03.75.Hh,03.75.Ss 

Recent experimental advances with trapped Fermi 
atoms enable one to reach considerably lower temper- 
atures than obtained previously, as well as to vary 
the effective attraction between Fermi atoms via Fano- 
Feshbach resonances^, ft then becomes possible to reach 
conditions where Cooper pairs of Fermi atoms form in 
weak coupling below the superfiuid critical tempera- 
ture Tc, and composite bosons form and Bose-Einstein 
condense in strong coupling, ft thus appears relevant 
to formulate a theory of the BCS-BEC crossover for 
trapped Fermi atoms for all temperatures in the broken- 
symmetry phase below Tc, connecting it with continu- 
ity to the results for the normal phase above Tc- This 
is the main task of the present paper, where a unified 
theoretical framework is set up for all temperatures and 
couplings. 

Limited theoretical results are so far available for 
the BCS-BEC crossover in a trap. A previous study 
of the density profiles over the whole crossover has 
dealt with the zero-temperature case within a mean- 
field approach^. Finite temperatures below Tc have been 
considered within mean field for a single coupling value 
in the weak-to-intermediate region"^. Fluctuations over 
and above mean field have been included over the whole 
crossover for temperatures above Tc^'^ 

In the present paper, we provide a systematic study 
of the whole BCS-BEC crossover in a trap by includ- 
ing fluctuations beyond mean field, for all temperatures 
below and above Tc, and up to the pair-breaking tem- 
perature T*. Below the critical temperature, our theory 
recovers the BCS results in weak coupling and the Bo- 
goliubov description for the composite bosons in strong 
coupling, and provides an interpolation scheme in the in- 
termediate (crossover) region where no small parameter 
exists for controlling the approximations of many-body 
theory. 

Study of the BCS-BEC crossover started with 
the pioneering work by Eagles for low-carrier doped 
superconductors^. A systematic approach to the prob- 
lem was later given by Leggett'', who showed that a 
smooth crossover from a BCS ground state of overlap- 
ping Cooper pairs to a condensate of composite bosons 
occurs as the strength of the fermionic attraction is in- 



creased. This study was later extended to finite tem- 
peratures above Tc by Nozieres and Schmitt-Rink with 
the use of diagrammatic methods* . Extension of this ap- 
proach to trapped fermions has relied so far mostly on a 
local Thomas-Fermi (TF) approximation'^''*'^. This local 
approximation is also adopted in the present paper. 

Our main results for the BCS-BEC crossover in a trap 
at finite temperature are the following: 

(i) We find that the critical temperature Tc increases 
monotonically from weak to strong coupling, reaching 
eventually the value of the Bose-Einstein temperature for 
the composite bosons in the trap. Correspondingly, no 
maximum is found in the intermediate-coupling region 
for the trapped case. The presence of this maximum was 
instead found with the same diagrammatic theory formu- 
lated for the homogeneous case®. 

(ii) We find that in the intermediate-to-strong coupling 
region the density profiles show a characteristic secondary 
peak located away from the trap center, at temperatures 
below but close to Tc. The occurrence of this peak is due 
to the combined presence of condensed and noncondensed 
composite bosons. We find that this peak survives up to 
couplings near the crossover region, such that the resid- 
ual interaction between the composite bosons is strong 
enough for the peak to be well pronounced. In this way, 
this peak could be experimentally accessible, providing 
one with a characteristic feature of the superfiuid state. 

(iii) We find that the "pairing fiuctuation" region be- 
tween Tc and T*, where precursor pairing effects should 
occur, is considerably reduced in the trap with respect 
to the homogeneous case. Pseudogap phenomena are 
thus expected to be very much reduced for trapped Fermi 
atoms, with respect to what occurs for high-temperature 
superconductors^ . 

The system we consider is a gas of Fermi atoms con- 
fined in a trap by a harmonic spherical potential V{r) 
(where r measures the distance from the trap center). 
The Fermi atoms equally populate two spin (hyper- 
fine) states and are mutually interacting via a point- 
contact (s-wave) attraction. This attraction is suit- 
ably regularized via the scattering length ap of the as- 
sociated (fermionic) two-body problem. The coupling 
strength is then identified with the dimensionless pa- 
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ramctcr (fc^a^)"'^, where the Fermi wave vector fc^? is 
related to the Fermi energy Ep = {BNy^^LO for non- 
interacting fermions in the trap^" by Ep = kp/{2m). 
Here, N is the total number of Fermi atoms, uj the 
trap frequency, and m the fermion mass (we set h = 1 
throughout). In principle, {kFaF)~^ w — oo corresponds 
to the (extreme) weak-coupling and {kpap)~^ « + oo 
to the (extreme) strong-coupling limit. In practice, the 
crossover between these limits occurs in the limited in- 
terval -1 < [kpap)-'^ < +1." 

The many-body diagrammatic structure for the ho- 
mogeneous case is considerably simplified by the use of 
the above regularization-'^^. In particular, in the broken- 
symmetry phase below Tc a diagrammatic theory for the 
BCS-BEC crossover can be set up^^ in the spirit of the 
t- matrix approximation'^^. This theory includes fluctua- 
tion corrections to the BCS results in weak coupling and 
describes the composite bosons in strong coupling by the 
Bogoliubov thcory-'^''. In the present paper, we extend 
this approach to the trapped case, by adopting a local 
TF approximation to take into account the trapping po- 
tential. This local approximation is implemented by re- 
placing the chemical potential /x, whenever it occurs in 
the single-particle self-energy and Green's functions, by 
the local expression /i(r) = /i — V{r). At the same time, 
the order parameter A is replaced by a local function 
A(r) to be determined consistently. 

Quite generally, in the BCS-BEC crossover approach 
the chemical potential is strongly rcnormalized when 
passing from the weak- to the strong-coupling limit. In 
our case, the coupled equations for the chemical potential 
and the local order parameter A(r) are: 

^(-) = -^/^ [iEft*..;MW.AW) 

m A(r) 



AT = 877 



l^e-^°^Gn(k,a;«;M(r),A(r)) 



(1) 



(2) 



where k is a wave vector and lUs = (2s -|- 1)tt/P 
(s integer) a fermionic Matsubara frequency {(3 being 
the inverse temperature). In the above expressions, 
5i2(k, tJs; /i(r), A(r)) is obtained from the BCS anoma- 
lous single-particle Green's function A/(£^(k)^ -|- w^) 
(with E{k) = yj(k)2TA2 and ^(k) = kV(2m) - 
fi), by replacing /z with /z(r) and A with A(r); 
the dressed normal single-particle Green's function 
Gii{\i,LOs; n{r), A{r)) is obtained by a similar replace- 
ment made on the function G'ii(k, w^), which contains 
fluctuation corrections beyond mean field. 

The quantity Gii(k, Wg) results from the solution of 
the 2x2 Dyson's equation in Nambu notation, with self- 



energy 

Ell ( k,Ws) = -E22(k, -Wg) 



'^^ ^ 5^ rii(q, 0^)^11 (q-k,f2,-a;«) (3) 



(27r)3 13 

Y,i2 ( h,ujs) = i;2i(k,u;,,) = -A 



(4) 



where = 27ri///3 [v integer) is a bosonic Matsubara 
frequency. Here, Qnik) = -(^(k) + iujs)/{Eik)^ + oi^) 
is the BCS normal single-particle Green's fimction and 
rii(g) = Xii(-«)/[Xii(9)Xii(-a) - xMqf] is the nor- 
mal pair propagator, with 

-Xn(9)=/(^ ^E^n(fc + g)an(-fc) 



(5) 



2E(k) 



Xi2(g) = /^;^E^i2(fc + g)a2i(-fc) 



(6) 



and the four- vector notation k = (k, ojs) and q ~ (q, fi^). 

Neglecting the diagonal elements (3) of the self-energy 
results in the BCS (mean-fleld) approximation. When 
extrapolated toward strong coupling, this approximation 
accounts for the formation of bound-fermion pairs upon 
lowering the temperature below T*. Inclusion of the di- 
agonal elements (3) of the self-energy is required to de- 
scribe condensation of these pairs at the lower tempera- 
ture Tc. In strong coupling, the normal pair propagator 
Fii (together with its anomalous coimterpart) reduce to 
the propagators for composite bosons within the Bogoli- 
ubov approximation. Above Tc, the diagonal elements 
(3) correspond to the t-matrix approximation in the nor- 
mal phase. 

Figure 1 compares the temperature vs coupling phase 
diagram for the trapped (t) and homogeneous (h) case, 
where Tc and T* arc identified (and normalized to the 
respective Fermi temperature Tp for the two cases) . The 
temperatures Tc and T* are obtained by solving the cou- 
pled equations (1) and (2) when A(r) = 0, with and 
without inclusion of the diagonal self-energy (3), in the 
order. 

Note that T|f increases monotonically from weak 
to strong coupling, approaching the value Tbe = 
0.9Auj{N/2y^^ of the Bose-Einstein temperature for the 
composite bosons in the trap'° (with the same trap fre- 
quency for fermions and composite bosons). No max- 
imum for Tc is thus found in the intermediate-coupling 
region for the trapped case, contrary to the homogeneous 
case where a maximum occurs at (kpap)^^ ^ 0.35. This 
behavior is consistent with the fact that, for a dilute Bose 
gas, interaction effects lead to a positive (negative) shift 
of the critical temperature in the homogeneous (trapped) 
case'^. Together with the vanishing of Tc in weak cou- 
pling, this implies that (at least) one maximum must be 
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present for the homogeneous case, while the presence of 
a maximum is not required for the trapped case. It is 
encouraging that our approximate theory leads to curves 
for the critical temperature in line with these general ex- 
pectations. 

Note further that the "pairing fluctuation" region of 
the phase diagram, delimited in each case by the curves 

T* and Tr, 



the density profile shrinks considerably and n{r = 0) 
increases correspondingly, as expected for a system of 
weakly-interacting (composite) bosons^'^. 
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FIG. 1. Temperature vs coupling phase diagram for the 
trapped (full lines) and homogeneous (dashed lines) CEise, with 

the critical temperature Tc and the pair-breaking temperature 
T* shown for the two cases. Each temperature is normalized 
to the respective Fermi temperature Tp. With this normal- 
ization the phase diagram is valid also for anisotropic traps. 

is considerably reduced in the trap with respect to the 
homogeneous case. In the strong-coupling limit, the re- 
duction of the pairing fluctuation region stems from the 
difference in the density of states D{c) at energy e for 
noninteracting particles, when passing from the homoge- 
neous {Dh{e) oc e-^/^) to the trapped {Dt{e) cx e^) case. 
This difference is, in fact, known to account for the larger 
value of {TBE/TF)t with respect to {TBE/TF)h^°. By a 
similar token, it can be shown that the same difference 
in the density of states accounts for the smaller value of 
{T*/TF)t with respect to {T* /TF)h- 

Figure 2 shows the density profiles n{r)/N (such that 
A-K fdrr'^n{r)/N = 1) vs r/Rp for three characteristic 
couplings, from T = to T* (where Rf = \/2Ef/ {muj'^) 
is the TF radius for noninteracting fermions in the trap). 
On the weak-couphng side of the crossover {{kFaF)~^ = 
— 1.0), n(r) depends mildly on temperature. In weak cou- 
pling, no information can thus be extracted from n(r) 
about when the superfluid phase is entered. For interme- 
diate couplings ((fci?aF)~^ = 0.0), n{r) starts to depend 
sensibly on temperature. Eventually in strong coupling 
{{kFaF)~^ = 1-0), n(r) shows a marked temperature de- 
pendence. In this case, the broad density profile at T* 
corresponds to a system of noninteracting fermions at 
the same temperature. Upon lowering the temperature. 
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FIG. 2. Density profile n{r)/N vs r/Rp for three coupling 
values and several temperatures from T = to T* . 

Note, in addition, the presence of a secondary peak 
away from the trap center, that shows up at temperatures 
below but close to Tc (this peak is most evident at about 
0.7Tc). The presence of this peak stems from the abil- 
ity of our theory to recover the Bogoliubov approxima- 
tion for the composite bosons. The presence of this peak 
was, in fact, predicted for a weakly-interacting trapped 
Bose gas^^. Our results show that this peak appears not 
only in the strong-coupling limit (which corresponds to 
weakly-interacting bosons), but also in the crossover re- 
gion. In this region, the residual interaction between the 
composite bosons is sufficiently strong for the peak to be 
well pronounced over the background, contrary to what 
occurs for weakly-interacting bosons. In this way, the 
presence of the secondary peak in n(r) below Tc could 
be subject to experimental testing, providing one with a 
characteristic signature of the superfluid state. 
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It is interesting to separate the total density n(r) for 
a given coupling into three components, namely, npir) 
for unbound fermions, no(r) for condensed pairs, and 
n'(r) for noncondcnscd pairs. These components are ob- 
tained from expressions similar to (2), with Gn therein 
replaced, respectively, by the Green's function for non- 
interacting fermions and by the differences Qn — Qo and 
Gil — Gil- By this procedure we project out from the 
many-body state its fermionic and bosonic character, not 
only in the extreme weak- and strong-coupling regimes 
where these components have independent physical real- 
ity, but also in the intermediate-coupling region of inter- 
est. 
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FIG. 3. Partial density profiles vs r/Rp for (kpaF)^^ ~ 0.0 
and T = 0.8Tc: nF{r)/N (full line), m{r)/N (dashed line), 
and n{r)/N (dotted line). 

The three components are plotted in Fig. 3 for the cou- 
pling value {kFaF)~^ = 0.0 and the temperature 0.8Tc. 
A layered structure results for these densities, from an 
inner core of Bose superfluid to an outer layer of normal 
Fermi liquid, in agreement with recent arguments for a 
dilute atomic- molecular Fermi cloud^^. Note that, for 
this intermediate coupling, the component n'(r) obtained 
by including fluctuations dominates over the condensate 
component no(r) that results from mean field. This im- 
plies that, already in the intermediate-coupling region, it 
is not possible to rely only on a mean-field calculation^ 
to obtain density profiles for trapped fermions at finite 
temperature. 

In conclusions, with a single theory that includes 
fluctuations beyond mean field, both below and above 
the critical temperature, we have studied the BCS-BEC 
crossover for a system of trapped Fermi atoms at finite 
temperature. Novel features, peculiar to the trapped 
case, have been contrasted with the results for the ho- 
mogeneous case. 
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